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We consider a class of single-particle one-dimensional stochastic equations which include external 
field, additive and multiplicative noises. We use a parameter 9 £ [0, 1] which enables the unification 
of the traditional Ito and Stratonovich approaches, now recovered respectively as the 6' = and 6 — 
1/2 particular cases to derive the associated Fokker-Planck equation (FPE). These FPE is a linear 
one, and its stationary state is given by a g-Gaussian distribution with q — ^!|!2m(i-9) ^ ^' where r > 
characterizes the strength of the confining external field, and M > is the (normalized) amplitude 
of the multiplicative noise. We also calculate the standard kurtosis ki and the g-generalized kurtosis 
ACq (i.e., the standard kurtosis but using the escort distribution instead of the direct one). Through 
these two quantities we numerically follow the time evolution of the distributions. Finally, we 
exhibit how these quantities can be used as convenient calibrations for determining the index q from 
numerical data obtained through experiments, observations or numerical computations. 

PACS numbers: 02.50.-r, 05.20.-y, 05.40.-a, 05.90.-|-m 



I. INTRODUCTION 

The random walk is the simplest model of diffusive 
processes in physics. If there is no wind, nor any other 
source of symmetry-breaking, a drunk has probability 
1/2 to take a step to the right, and the same probability 
to take a step to the left, at each instant. Underneath 
this model there is a formal mathematical theory known 
as stochastic calculus, which in some sense can be inter- 
preted as a extension of the standard differential calculus 
taught in undergraduate courses. The stochastic as well 
as the standard calculus are based on the definition of 
the integral. Let us define the stochastic integral by 



I[G{t)] = / dW{t')G{t') 



(1) 



where G{t) is a left-continuous function (i.e., a function 
which is continuous from the left at all the points where 
it is defined) and W{t) is a Wiener process [H, Ql- As 
in the Riemann integral definition, the formal stochastic 
integral (also known as Riemann- Stieltjes integral) is a 
infinite discrete sum of very small intervals {dW{t')) of a 
stochastic function. When we perform this sum we must 
make a choice, more precisely the function G{t') has to 
be evaluated in some point inside each interval. This 
choice defines what kind of stochastic calculus will be 
performed henceforth. The two most famous procedures 
are Ito calculus and Stratonovich calculus. In the first, 
Ito calculus, the function is evaluated at the beginning 
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of the intervals: 

n 

Ii[G{t)] = ms - lim y G,_iAWi , (2) 

i=l 

where (ms-lim) stands for mean squared limit, i.e. a 
second moment convergence; n is the number of subin- 
tervals in which we divide the interval [to;^]i ^W^j = 
Wi — Wi-i. In the Stratonovich calculus we take the 
arithmetic average between the integrate function values 
at the beginning and at the end of the intervals, as follows 



Is[G{t)] = ms-lim y 



-AW,, 



(3) 



i=i 
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An unified form can be proposed [3|-|6| for the stochas- 
tic integral, namely 

n 

Ie[G{t)] = ms - lim V [6G, + (1 - 0)G,_i] AW, , (4) 

71 — Vr^'j < ' 



where < 6* < 1. Notice that the two traditional pro- 
cedures can be recovered easily. Indeed, if 6* = we 
obtain the Ito approach, and if = 1/2 we obtain the 
Stratonovich approach. Moreover, li 9 = 1 we recover 
the so called backward-Ito Q stochastic approach, that 
also known as isothermal convention 

It is possible to go one step further and generalize 
(HI). Indeed, we can assume that the values oi 9, at 
each interval dW{t'), are given by an arbitrary distri- 
bution p{9) (Jg d9 p{9) = 1). In particular, if this dis- 
tribution is p{9) ~ 5{9 — 9a), we recover the three cases 
mentioned above for suitable values of 9q, namely = 
(Ito), 9o = 1/2 (Stratonovich) and 9o = 1 (backward- 
Ito). An interesting remark arises when p{9) is constant. 
In this case (9) = 1/2 which coincides precisely with the 



2 



average value corresponding to the (frequently adopted) 
Stratonovich approach. In some sense, this is what seems 
to happen in most experiments, where the act of measur- 
ing is not instantaneous but in a time window. Let us 



J 



finally mention that, in some sense, the only hypothesis 
strictly consistent with causality 8] is p[0) — 5{9). 

An instructive example of the aplication of definition 
dl]) consists in the choice G{t) = W{t), so 



Ig[W{t)]= / dW{t')W{t') ^ ms-lim V [6*1^4 + (l-6')W^^_i]AW^i 

tn A — 1 



(5) 



r 



Again, when 6* = (0 = 1/2) we obtain the standard re- 
sult for the Ito (Stratonovich) approach. For more details 
about the derivation of the last equation, see Appendix. 

In addition to the ingredients of stochastic calculus 
that we have mentioned above, some other concepts will 
be necessary in the present paper. In particular, escort 
mean values [3, 111] are convenient theoretical tools for 
describing basic features of some probability densities, 
mainly those densities which decay as power laws that 
naturally appear in the study of complex systems dy- 
namics, such as those obeying nonextensive statistical 
mechanics 16]. This theory generalizes the Boltzmann- 
Gibbs (BG) statistical mechanics, and is governed by an 
entropic index q, which equals unity for the BG case. 
The characterization of a probability density by its set of 
escort mean values, if all of them converge, is a natural 
extension of the well-known characterization of a distri- 
bution in terms of its standard moments, which corre- 
sponds to q = 1. The g-generalized theory has been 
applied to calculate many features of several complex 
systems, such as: (i) The velocity distribution of cells 
of Hydra viridissima, that follows & q — 3/2 PDF [l7| : 

(ii) The velocity distribution of (cells of) Dictyostelium 
discoideum, that follows a q ~ 5/3 PDF in the vegeta- 
tive state and a. q = 2 PDF in the starved state 18|: 

(iii) The velocity distribution in defect turbulence 19|; 

(iv) The velocity distribution of cold atoms in a dissipa- 
tive optical lattic e |22| : (v) Velocity distribution during 
silo drainage |23, |2l[; (vi) The velocity distribution in a 
driven-dissipative 2D dusty plasma, with q — 1.08 ±0.01 
and q = 1.05 ± 0.01 at temperatures of 30000 if and 
61000 K respectively [11]; (vii) The spatial (Monte Carlo) 
distributions of a trapped ^^^J5a+ ion cooled by various 
classical buffer gases at 300 K [23| ; (viii) The distribu- 
tions of price returns at the stock exchange [25l - l27i] ; (ix) 
The distributions of returns of magnetic field fluctuations 
in the solar wind plasma as observed in data from Voy- 
ager 1 ^] and from Voyager 2 [29] ; (x) The distributions 
of returns of the avalanche sizes in the Ehrenfest's dog- 
flea model [131; (xi) The distributions of returns of the 
avalanche sizes in the self-organized critical Olami-Feder- 



Christensen model, as well as in real earthquakes 
(xii) The distributions of angles in the HMF model 



31| 



(xiii) The distribution of stellar rotational velocities in 
the Pleiades [lH ; (xiv) The distribution of transverse mo- 
menta in high-energy proton-proton collisions [sij; (xv) 
In the ation of spin glasses [3a] . 

In the present work we use Eq. ([5]) to write unified 
Langevin and Fokker-Planck equations. In section |lTl af- 
ter deriving these equations, we will show that the steady 
solutions of the Fokker-Planck equation (FPE) are q- 
Gaussians (see later on for their precise definition) whose 
entropic parameter q depends on the noise and drift am- 
plitudes, and also on the choice of a stochastic approach 
represented by the parameter 9. In section Hill we intro- 
duce a generalized kurtosis based on [Tij to characterize 
the temporal evolution of the distribution. After that, 
we integrate numerically the FPE with an initial distri- 
bution different from its asymptotic form. In particular, 
we consider as initial distributions g-Gaussians charac- 
terized by an index qi [qi ^ q). Our numerical results 
show how the convergence towards the attractor behaves 
as a function of the parameter 9. Finally, we compare 
the standard kurtosis {q = 1) with the one calculated 
using escort mean values, namely g-kurtosis. Our re- 
sults show that the standard kurtosis has a divergence 
at q — 7/5 while the g-kurtosis has no divergence in the 
range — 1 < g < 3. In addition to that, we show that the 
standard kurtosis and the g-kurtosis are monotonic func- 
tions of the entropic parameter q, which suggests that 
they could be used as calibration curves to determine, 
from numerical data, the most appropriate value of q. 
Let us now describe the consequences of the unification 
on the Langevin and Fokker-Planck descriptions. 



II. FOKKER-PLANCK EQUATION AND ITS 
SOLUTIONS 

A stochastic differential equation (SDE) is not com- 
pletely defined by itself. If the Langevin equation has 
multiplicative noise, we must choose of what approach, 
Ito or Stratonovich, will be used to integrate it. This is 
the well-known Ito-Stratonovich dilemma [11, [2, [13 and 
is ultimately solved by taking into account the specific 
features of the system under investigation. For instance. 
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if the noise has a finite correlation time Tc (even if the 
Umit Tc — !> is used to derive the SDE) or the noise comes 
from external sources [s^ , the Stratonovich choice is the 
adequate one. Instead, the Ito formalism is the correct 
choice if T>: is strictly zero or the noise comes from internal 
sources [H, [l^l ■ A very interesting and enlightening dis- 
cussion about the Stratonovich-Ito dilemma can be found 
in This controversy between the two approaches 

appears in the form of the so called noise-induced drift 
which is an effect of the state dependence of the noise 
strength. 

Based on the proposal (U) we will write an unified 
Fokker-Planck equation that contains, as special cases, 
the Ito {9 = 0) and Stratonovich {9 = 1/2) forms. To do 
it, let us consider the quite general SDE for a stochastic 
variable u{t): 



du{t) ^a[u{t),t\ dt + b[u{t),t\dW{t) 



(6) 



where a[u{t), t] is a deterministic external force given, for 
instance, by a potential, h[u{t),t] is the state-dependent 
noise amplitude, and dW{t) = ^{t)dt is a Wiener incre- 
ment (^(i) is a Gaussian, zero mean, white noise). We 
will use a shorter notation for the dynamic variable u(t) 
omitting its time dependence, so from now on u = u{t), 
a[u{t),t] = a{u,t) and b[u{t),t] = b{u,t), unless it is in- 



dispensable for the clarity of the text. 

We are interested in the dynamics of the distribu- 
tion function of u, so once we have the Langevin equa- 
tion we may obtain a FPE for the probability den- 
sity P{u,t) by the Kramers-Moyal expansion Q dtP = 
E„>i(-5«)"[^'"-'^]> where the coefficients D^") are 
given by 



(a;,t) = -7 lim \u{t)=x ■ (7) 

n! c-*o e 



To calculate these coefficients, we need to write the 
Langevin equation in the integral form 

pt+t rt+c 

u{t+e)~x^ dt'a[u{t'),t']+ dW{t') b[u{t'),t'] , 
Jt Jt 

(8) 

and assume that a{u,t) and b(u,t) can be expanded as 

a[u{t + e),t + f] = a[x,t + f] + a'[x,t + f]e+... (9) 
b[u{t + e),t + f] = b[x,t + f] + b'[x,t + f]e+ {10) 

where a' and 6' means differentiation with respect to u. 
Putting dHl) and Uni), up to first order on e 0|, into (0) 
and iterating the result, we get 



{u{t + e)-x) = ea{u,t + e)+b\u,t + ()b{u,t + e)l^j^ W{t')dW{t')^ 

= e[a{u,t + e) + 9b'{u,t + e)b{u,t + e)\ . 

I 



(11) 
(12) 



If the noise in the Langevin equation is (5-correlated, after 
repeated iterations it is possible to show that, at e-order, 
the terms of the Kramers-Moyal expansion with n > 3 
vanish. Using this arguments we can calculate the first 
two coefficients of that expansion, namely 

Kg{u.t) = D'^^Xu.t) = a{u,t) + 9b'{u,t)b{u,t){li) 
D{u,t) = D^^Hu,t) = b^{u,t) . (14) 

This procedure leads to the following Fokker-Planck 
equation: 

dP 5 r_ , , „ Id 



dt 



du 



Ke{u,t)P^ 



2du 



[Diu,t)P] 



(15) 



We recover, for 9 — and 9 — 1/2, respectively, the 
Ito and Stratonovich forms currently found in the lit- 
erature. As expected, the difference between the two 
forms appears in the drift term Kg, not in the diffusion 
term. If we deal with a FPE in the Stratonovich form 
we should include the above mentioned noise-induced 
drift term ^du[D{u)], which is unnecessary in the Ito 
form. We can write these FPE as a continuity-like equa- 
tion dtP = —duj{u), if we define j{u) = Kg{u)P — 



^du[D{u)P] as a probability current. In this case, the 
stationary solution of the FPE can be obtained ele- 
gantly from zero flux boundary conditions j(±oo) = 
j{u) = 0. It is very interesting to notice that this prob- 
ability current can be rewritten as j{u) — Kq{u)P — 
\[D{u)f'^du{[D{u)Y-'^'^P], which allow us to rewrite 
the FPE as: 

(16) 

Notice that here, in contrast with Eq. (fTSj). 9 only ap- 
pears in the diffusion term, and not in the drift one. The 
demonstration of the equivalence between Eqs. ^5]) and 
(1161) is straightforward and can be found in the Appendix 
llYBl 

Let us now consider a family of models represented by 
Langevin equations of the type: 

f{u)+9{u)m+v{t) , (17) 
where ^(i) and 7]{t) are uncorrelated zero- mean Gaussian 
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white noises with autocorrelation function given by: 

mm) = 2A/<5(t - t') , {vitUt')) - 2A6{t - t') , 

(18) 

where M > and A > 0. The Eq. dTT]) can be rewritten 
as: 

f{u)+g{u)at) , (19) 

where the additive and multiphcative noise terms were 
replaced by an effective multiplicative noise given by 
g{u) = \J (Af[(jr(u)]2 + v4)/C and C(0 being a zero-mean 
Gaussian white noise with autocorrelation function given 
by (C(t)C(i')) = 2C(5(i - t') (C > 0). A possible demon- 
stration of this relation is given in the Appendix. No- 



tice that Eq. (HH) is of the form of Eq. ([6]) with 
dWit) = C{t)dt. 

As we want to analyze the influences of the additive 
and multiplicative noises separately, we will use the Eq. 
(IT71) instead of (HH) henceforth. 

Following the same procedure adopted to calculate the 
coefficients of the FPE associated with Eq. (jG)), we can 
easily calculate the Kramers-Moyal coefficients for (|17|) . 
They are: 

Keiu) = f{u) + 29Mg{u)g'{u) (20) 
D{u) = A + M[g{u)f (21) 

These results lead us the following FPE: 



dP{u,t) 
Ft 



du 



{[f{u) + 20Mg{u)g'{u)] P{u, t)} + M-^ { [giu)f P{u, t)] 



A—P{u,t) 



(22) 



sons suggesting that, in the extremization of Sq , it is 



which is basically Eq. (jlSp with an additional term due propriate constraints. There are many theoretical rea- 
to additive noise. 

For f{u) derived from a potential-like function V{u) = 
(t/2) [(;(u)]^, the stationary solution has the form 
P(u, oo) cx eq ^^^'^'^ Q where the g-exponential function 
is defined as follows: = [1 + (1 — q)x\^_l^-^^''\ with 



[z\+ = z if z > 0, and zero otherwise; ef — e^. If we 
consider the simple case g{u) oc u, then the stationary 
state distribution is a g-Gaussian, i.e., P{u, oo) cx e~^" , 
where q and (3 are given by: 



convenient to express the constraints that are being im- 
posed in a escort mean value form. Escort mean values 
(or g-moments) are useful tools to analyze power law dis- 
tribution that frequently appear in the study of complex 
systems. They are defined as 



/3 = -7]- (23) 

^ r + 2Af(2-^) 
^ r + 2M(l-6l) ^ ' 

We can verify that these results unify those obtained in 
0, more specifically q = ^'^2M ^'^^ '^^^'^ = 0)j 

and q = '^^^ for the Stratonovich case {9 = 1/2). 

Let us remind at this point that the g-Gaussian form 
precisely is the one which, under appropriate constraints, 
extremizes the entropy 

1 - / dx[p{x)]i 



Sq[p\ = ks- 



q-1 



(25) 



where Si = Sbg — J dxp{x) In p{x). For further 
connections between the structure of Fokker-Planck-like 
equations and entropy see [H, [13] ■ 



III. GENERALIZED MOMENTS AND 
KURTOSIS 

As already mentioned, the g-Gaussian is the distribu- 
tion form which extremizes the entropy (|25p under ap- 



-\-oo 



Aix)fqix)dx, 



where the escort probability density is given by 



j:::[m]'^dx 



(26) 



(27) 



The characterization of a probability density in terms 
of its escort mean values is a natural extension of the 
well-known characterization of a distribution in terms of 
its standard moments if all of them are finite. The phys- 
ical interpretation of 5-generalized mean values demands 
an explanation. The quantities that are physically im- 
portant in order to characterize a distribution are, for in- 
stance, the most probable value, the range of the typical 
values, the width of the distribution, its asymmetry, and 
so on. Such information is conveniently contained in the 
set of successive mean values of the distribution, as long 
as they are finite. What can be done whenever all mo- 
ments above a given one diverge? For example, if we are 
dealing with the Cauchy-Lorentz distribution, how can 
we characterize its width (obviously finite for any given 
such distribution) ? Certainly not through its second mo- 
ment, since it diverges! By appropriately choosing the 
value of q (see below), its width can be characterized by 
its g- variance (i.e., its variance with the escort distribu- 
tion), which will also be finite. 
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In [37|, the correct set of all escort mean values is 
shown, together with the set of all associated normalizing 
quantities, that characterize a given probability density 
fix), even if it decays as slowly as a power-law. Based 
on a generalization of the Fourier transform, namely the 
q-Fourier transform, defined by jss! ] 



/ + 00 
-OO 



dx. 



(<?>!), (28) 



it is possible to expand the g-characteristic function and 
obtain the correct exponents values necessary to perform 
the calculations of the g-moments. The family of escort 
mean values which arises from this procedure are given 
by 



where 



n^^''[f{x)v-dx 

J+^[f{x)]'^^dx 



l + n(g-l) 



(29) 



(30) 



The kurtosis, usually defined in terms of the ratio be- 
tween the fourth (n — 4) and three times the squared sec- 
ond moment [n = 2) of a given distribution, is regarded 
as a measure of how different a given distribution is from 
a Gaussian. Therefore the g-kurtosis, defined in terms of 
g-moments, is the analogous measure for g-Gaussians. Its 
mathematical forms are defined, respectively, as follows: 



j:^^x^[mr-'dx 



•{!^^x^fix)f^-'dxy 



(31) 



Notice that the standard kurtosis can be regarded as a 
particular case (g = 1) of Eq. (|3T|) . Notice also that 
the number 3 in the denominator of Eq. pip plays no 
special role, and can equally well be replaced by the q- 
dependent value exactly corresponding to g-Gaussians, 
instead of that for Gaussians. 

To calculate the correct value of g we shall now con- 
sider a typical situation arising in complex systems, such 
as those commonly addressed within g-thermostatistics. 
Usually these systems asymptotically behave as power- 
law probability densities, i.e.. 



fix) 



!| — >■ oo; 7 > 1) . 



(32) 



It is clear that, if fix) is not defined on a bounded in- 
terval, the standard linear moments (a;") may diverge 
above some value of n. Therefore the standard proce- 
dure to characterize the probability density through all 
its moments cannot be implemented. However, in ,37 ] it 
is shown how the escort mean values can overcome such 
difficulties. In this work, the relation is established be- 
tween the exponent 7 of the power-law distribution, and 
the escort mean values parameter g. It is given by 



1 



1 



(33) 



So, once we know the power-law exponent 7, we will have 
the correct value of g, hence of the values g„ which enable 
the calculation of the escort mean values of the probabil- 
ity density. 

For the special case where fix) — eq '^ (a Q- 
Gaussian), we have: 

eg''"' - (|x|-^oo), (34) 

hence 7 = 2/((5 — 1), for Q > 1. In this case, the relation 
became: 



9-1 = 



1 



(35) 



A Q-Gaussian is normalizable for Q < 3; for Q < 1 it 
has a compact support, and an unbounded support for 
1 < Q < 3. Its second and fourth moments diverge for 
5/3 < Q < 3 and 7/5 < Q < 3, respectively. But their 
second and fourth Q-moments are finite for Q < 3. These 
two Q-moments can be analytically calculated [1^, l37j 
and written in terms of Gamma functions. 



)lQ-\ 



2(Q-1)^^(q^-|) 

2{i-Q)r{^+^) 
4(Q-i)=r(^+i) 

3/j-^r(^+i) 
4(i-Q)'/^(t^"+|T 



for 1< Q < 3 

for Q = l (36) 

for < Q < 1 

for 1 < g < 3 
for g = 1 (37) 

for < g < 1 . 



We intend to compare the temporal evolution of the 
distributions given by the FPE of the previous section 
to its stationary states, which we already know to be Q- 
Gaussians. To do it, it is convenient to normalize the Q- 
kurtosis by its stationary value, k,q{co) = limf^oo Kqit), 
which is given by Eqs. (PSI) and ((57|) . 



i^gioo) = < 1 



IT 



^^(Tjg + l)^(T% + §) 



for 1< g < 3 
for Q = 1 

for < g < 1 



(38) 

In Fig. [1] we can see the temporal evolution of the nor- 
malized g-kurtosis obtained by the numerical integration 
of the FPE. In the left column, we show the temporal 
evolution for the three stochastic procedures above men- 
tioned, Ito i9 ~ 0), Stratonovich (^^ = 1/2) and backward 
Ito i9 — 1). In the right column we show the same tempo- 
ral evolution of the normalized g-kurtosis but now with 
6 = and different values of Qi. As expected, the kurto- 
sis approaches monotonically its stationary value for all 
values of parameter Qi. 

In Fig. [5] we can compare the behavior of the standard 
kurtosis and the g-kurtosis given by Eq. (pS)) . As we 
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FIG. 1: On the left column, the temporal evolution of Q- 
kurtosis for Ito (0 = 0), Stratonovich {6 = 1/2) and backward 
Ito {9 — 1) stochastic procedures. The numerical results are 
normalized by its stationary value Kq = Kqipo). In (a) we 
have Q = I and Qi = 5/4, in (b) Q = 7/5 and Qi = 1.65 
and (c) Q = 1.666... and Qi = 1.91666... . As expected, 
for sufficient large times the numerical results approaches to 
unity. The right column shows the temporal evolution of the 
normalized Q-kurtosis with a fixed 9. In (a) we have (5 = 1 
and Qi = 1, 5/4 and 3/2; (b) shows the results for Q = 7/5 
and Qi = 1.15, 7/5 and 1.65; and (c) the results for Q = 5/3 
and Qi = 1.41666..., 5/3 and 1.91666... . 



can see, the Q-kurtosis diverges at Q = —1, but is finite 
for — 1 < Q < 3, while the standard kurtosis diverges for 
7/5 < Q < 3. As expected, the Q-kurtosis coincides with 
the standard one for Q = 1. Furthermore, the Q-kurtosis 
is a monotonically decreasing function of Q. Based on 
this feature, experimentalists could use it as a calibration 
curve that allows the determination of the proper value 
of the entropic index Q for a given system. 
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FIG. 2: The standard kurtosis ki and the Q-kurtosis kq as 
functions of Q. We verify that both have the same value 1 
for (5 = 1. For Q > 1, the use of kq is generically more 
convenient than using k\. Indeed, the latter one diverges at 
Q — 7/5 while kq remains finite up to the maximal value 
(5 = 3. In contrast, for Q < 1, ki is more convenient than 
KQ. Indeed, kq diverges at Q — —1, whereas ki remains finite 
down to (5 — > — oo. In the Q — oo limit, ki saturates at 
the value 3/5. 



plicit form for the Q-generalized kurtosis to study the 
convergence towards the Q-Gaussian stationary solution. 
The difference in the convergence for the same q-Gaussian 
with different 0's is due to the change in the amplitude of 
the multiplicative noise needed to compensate the change 
in 9 (see Eq. ([24])). This is the origin of the difference 
in the convergence velocity when we keep Qi's and Q's 
constant, and vary 6. Looking at the time evolution of 
the Q-kurtosis starting from different initial conditions, 
i.e. different Qi's, we see that the convergence is reached 
for a little less than two decades in all explored cases. 
Finally, we propose that, if a system is well described 
by the nonextensive statistical mechanics, and has a Q- 
Gaussian form for its probability density function, the 
Q-kurtosis can be used as a calibration curve to deter- 
mine, from data, the best value of the entropic index Q. 
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Appendix 

A. The integral Ie[W{t)] = //^ dW(f')W(t') 

Here we perform, in details, the calculation of the 
stochastic integral of the function G{t) = W{t) in terms 



7 



of discrete sums. 



The definition of the unified stochastic integral is 



Ie[G{t)] = ms - Km 5„ = ms - lini V [OG^ + (1 - 6i)G,_i] AWi . 

r>. — >-rxn n. — ^i-io * ^ 



(39) 



r 



If G(t) = W{t), we have 

n 

Sn=J2 + (1 - 0)W^,_iAW^,] . (40) 

1=1 

The terms WiAWi and Wi-iAWi can be written as 



under a summation. Therefore only the first and 
the last terms need to be taken into account in the sum- 
mation. Indeed, all the others mutually cancel, and we 
will finally have 



(43) 



2W,AW, = {W,f + {AW^f ~ {W, ~ AW,f (41) y^^ ^^^^ ^ W^jt) -WHt,) -Y^Ui^^W^f ^^^^ 



Wi 



Notice that in both equations we have terms like {Wf 



Taking this last two equations into account and the fact 
that ms -lim„^oo Yll=i{^^if' = t ^toi obtain 



Ie[W{t)] = ms - hm 5„ = - {[W{t)]'' - [Wito)? - (1 - 29){t - to)} 



(45) 



which is the same result presented in Eq. 



left hand side. 



B. How to get Eq. ([15]) from Eq. (fT6\i 

In section |ll] we claim that Eqs. (jTSj) and ([TB]) are 
equivalent. To prove this statement, we only need to 
show the following equality: 

where the dependence on the dynamical variable u was 
omitted. To do this, we perform the derivative on the 



D 



20 



So, 



{l-29)D-^'P^+D^-^'^ 
ou ou 



Ou 



OU ou 



(46) 



J 



KoP--[Dr-^{[D]''''P} 



If dD d , , 



(47) 



r 



where K^P + OduDP = Kq. 



C. The equivalence between ii — f{u) + giu)S,{t) + ri{t) 
and ii = f{u) +g(it)C(i) 



Suppose two different Langevin equations for the 
stochastic processes u{t), both with the same determin- 
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istic drift force f{u) but different noise sources, Let us calculate the variance corresponding to both equa- 

tions. From the first one we obtain 

u = f(u) + g{u)^{t) + r^it) (48) 
ii = f{u) + ~g{u)C{t) . (49) 

The mean- value solution of both equations is 

(u(i)) = fdt'fiu) . (50) 
I 

{[uit) - {uit))f) = dt' £ dt" {g[u{t')]i{t') + ry(t')} {9Ht")W) + 7?(t")}) (51) 

'dt' fdt" { g[zi(iO]5[«(t'0] (£2252!2+^2^£i£2l ^^^^ 

^2MS{t'-t") =2AS{t'-t") 



-'0 



t 

dt' 







+g[u{t')]{mf''))+9Hn]{my{t'))} (53) 
=0 =0 

' {2M{g[u{t')]f + 2Aj . (54) 



From the second one we obtain 



{[u{t) ~ {u{t))f) = J' dt' 1^ d^"ff[u(^0]g[z.(^")] (C(^OC(^")) 

2CS(t'-t") 

dt'2C{g[u{t')]f . 



(55) 
(56) 



If we want that both Langevin equations give us the same 
mean-squared value for the stochastic process u{t), we 
must match the integrands of Eqs. ([?T|) and This 
leads to the following relation among noise terms: 



This result shows that Eqs. like ((T7)) can be rewritten as 
in (jl9p . with an effective multiplicative noise term given 
byEq. (E 



M[g{u)]^ + A 
C 



(57) 
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